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LONG-TERM  GOALS 

(a)  To  provide  an  accurate,  efficient,  computational  model  (two-scale  approximation,  TSA)  for  the 
4-wave  interactions,  in  operational  wave  forecast  models,  suitable  for  global,  basin  and  coastal 
scale  applications,  and  able  to  transition  seamlessly  from  deep  to  shallow  water. 

(b)  Fully  test  TSA  with  respect  to  exact  codes  for  the  full  Boltzmann  integral  (FBI),  for  duration- 
limited,  fetch-limited  wave  growth,  turning  winds,  swell-windsea,  interactions,  etc. 

(c)  Numerically  investigate  and  clarify  the  basis  for  TSA,  its  limitations,  errors,  enhancements, 
improvements,  self-similarity  properties,  and  spectral  flux  properties. 

(d)  Implement  TSA  in  a  variety  of  modern  operational  wave  forecast  models,  e.g.  WAVEWATCH™ 
(WW3)  and  SWAN  for  extensive  tests  on  important,  realistic  wave  conditions. 

(e)  Derive,  adapt  and  implement  new  formulations  for  source  terms,  wind  input  Sin,  and  dissipation 
Sds,  from  recent  literature  and  the  NOPP  partnership,  with  TSA,  in  modern  wave  models,  for 
tests,  including  veering  or  accelerating  winds,  sea  and  swell  interactions,  and  real  storm  cases. 

OBJECTIVES 

For  this  reporting  period: 

1)  Enhance  the  formulation  of  TSA,  for  computational  efficiency,  and  also  to  consider  complicated 
situations,  such  as  turning  winds,  and  real  storms,  implemented  in  WAVEWATCHIII  (WW3). 

2)  Complete  tests  with  storm-generated  waves  for  TSA,  using  new  Sin  and  Sc/s  paramctcrizations,  in 
WW3,  for  real  storms,  as  well  as  hypothetical  fetch-  and  duration-limited  wave  growth  situations. 

3)  Compare  TSA  implemented  in  WW3,  including  tests  with  relatively  simple  parameterizations  for 
Sm  and  Sds  forms  for  wave  models.  Include  tests  for  SWAN  for  shallow  water  waves  cases. 

4)  Finalize  the  TSA  formulation  in  terms  of  computational  efficiency  and  accuracy. 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2Q 13  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2013  to  00-00-2013 

4.  TITLE  AND  SUBTITLE 

TSA  -  a  Two  Scale  Approximation  for  Wind-generated  Ocean  Surface 
Waves 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROTECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  North  Florida, School  of  Engineering,!  UNF 

Drive, Jacksonville, FL, 32224 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBIECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  OS 

unclassified  unclassified  unclassified  Report  (SAR) 

6 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


APPROACH 


We  will  clarify  the  limitations  and  properties  of  the  finalized  TSA  implemented  in  modem 
operational  wave  models,  with  consistent  source/sink  tenns,  Sin  and  S*. 

1)  Enhanced  versions  of  TSA.  We  will  complete  tests  and  derivations  of  new  optimized  versions  of 
TSA,  aiming  to  reduce  execution  times,  maintain  or  improve  accuracy,  and  address  areas  where 
TSA  has  been  shown  to  be  deficient,  compared  to  FBI  (full  Boltzmann  integral),  e.g.  accelerating  or 
veering  winds,  combined  sea/swell  conditions,  etc.  This  activity  must  address  any  TSA  deficiencies 
that  have  become  evident,  in  this  project. 

2)  TSA  tests  with  new  S,„  and  S*.  New  parameterizations  for  Sin  and  S*  and  modifications  to  TSA 
necessitate  the  need  to  redo  basic  tests  for  duration-limited,  fetch-limited,  turning  winds,  and 
additional  tests,  using  TSA  in  WW3,  as  are  always  needed  to  ensure  that  the  model  is  reliable. 

3)  TSA  model-model  comparisons.  We  will  compare  selected  cases  for  TSA  implemented  in  both 
WW3  and  SWAN  with  corresponding  model  estimates,  using  the  variety  of  accepted  formulations 
for  Sin  and  S'*,  including  new  forms  as  in  item  2),  and  older  parameterizations,  for  stonn  cases, 
comparing  model  results  with  observed  data. 

4)  Finalized  TSA  and  source  terms.  The  TSA  formulation  derived  in  this  project,  will  be  finalized, 
optimized  in  terms  of  run  time  and  efficiency,  and  implemented  in  recent  WW3  and  SWAN 
versions.  This  includes  additional  testing  and  optimization  of  TSA  for  the  new  fonnulations  of  Sin 
and  S*  to  get  acceptable  results  compared  to  observed  data.  Shallow  water  tests  are  included.  Tests 
for  reliable  modem  relatively  simple  operational  wave  models  are  completed  by  Resio,  separately. 

WORK  COMPLETED 

The  project  has  completed  the  objectives  that  were  originally  proposed. 

1 )  Enhanced  versions  of  TSA.  In  evolving  wind  and  wave  conditions,  the  TSA  formulation  works  well 
for  simple  fetch-limited  or  duration-limited  wave  growth.  However,  when  the  wind  direction 
continuously  changes,  the  TSA  fonnulation  needs  to  be  modified  to  allow  the  broad-scale  term  to 
respond  properly.  If  there  is  a  misalignment  of  TSA’s  broad-scale  tenn  with  respect  to  the  change  in 
wind  direction,  and  TSA  is  not  able  to  provide  a  reliable  representation  of  the  nonlinear  transfer  due 
to  wave-wave  interactions  in  this  changing  situation.  The  multiple-TSA  approach  allows 
generalization  of  the  broad-scale  tenn,  allowing  more  than  one  broad-scale  parametric  tenns, 
conesponding  to  multi  -peaked  spectra.  It  can  be  shown  that  this  modification  of  TSA  (denoted 
‘double  TSA’,  or  dTSA)  allows  the  TSA  approach  to  respond  as  the  wind  direction  veers,  and  is 
important  for  situations  which  result  with  swell-  windsea  interactions.  Thus,  model  accuracy  is 
improved,  for  waves  generated  by  storms  simulations.  For  additional  details  see  Perrie  et  al. 

(2013a;  2013b;  2014),  and  Figure  1  below. 

2)  TSA  tests  with  new  Si„  and  Sds  parameterizations .  We  have  completed  comparisons  with  modem 
source  terms  ( Ardhuin  et  al.,  2010),  so-called  ‘ST4’,  in  comparison  with  older  WAM  cycle  physics, 
from  WAM  model  (WAMD1,  1988),  and  alternate  formulations,  by  Tolman  and  Chalikov  (1996), 
and  others,  and  to  investigate  the  response  of  TSA  to  these  source  terms,  in  stonn  conditions  in  the 
North  Atlantic.  Without  any  parameter  tuning  for  ST4  (that  is  Sin  and  S*),  we  show  that  wave 
model  results  from  the  implementation  of  TSA  in  WW3  are  better  than  results  obtained  by  earlier 
source  term  formulations  for  Sin  and  S'*,  or  mns  using  DIA  -  the  discrete  interaction  approximation. 
See  Perrie  et  al.  (2013a;  b;  2014). 


3)  TSA  model-model  comparisons.  TSA  was  implemented  in  recent  versions  of  WW3  (4.10  and  4.18), 
and  SWAN,  and  tests  were  made  to  compare  model  skill,  using  modem  source  terms  of  Ardhuin  et 
al.  (2010),  denoted  ‘ST4’,  as  well  as  simpler  source  terms,  like  ST1  (WAMD1,  1988).  Tests  include 
simple  hypothetical  cases,  as  well  as  storms  such  as  hurricane  Juan  (2003),  and  mid-latitude 
nor’easters.  We  show  in  these  model  studies  that  TSA  works  well  with  modem  source  terms  like 
ST4,  and  that  model  skill  is  improved,  compared  to  older  versions  of  these  source  terms.  Results  are 
shown  in  Perrie  et  al.,  (2013a;  2013b;  2014). 

4)  Finalized  TSA  and  source  terms.  The  derived  TSA  formulation  was  optimized  in  terms  of  mn  time 
and  efficiency,  and  implemented  in  WW3  and  SWAN.  This  involves  an  optimized  selection  of  grids 
and  computation  loops  over  which  the  computation  is  carried  out.  The  SWAN  model  tests  use  an 
unstructured  grid  allowing  high  accuracy  for  coastal  areas  with  shallow  water.  Shallow  water  tests 
are  included  in  tests  with  unstructured  SWAN,  including  hurricane  Juan  (2003).  Results  are 
reported  by  Sun  et  al.,  (2014).  Additional  tests  using  a  reliable  modern  simple  operational  wave 
model  were  completed  as  part  of  a  separate  project  led  by  Resio. 

RESULTS 

The  meaningful  technical  results  achieved  in  this  fiscal  year  are: 

a)  We  have  developed  an  enhanced  version  of  TSA  which  is  a  generalization  of  the  standard  TSA 
formulation,  originated  by  Resio  and  Perrie  (2008),  and  Perrie  and  Resio  (2009)  in  the  sense  that  it 
is  able  to  not  only  simulate  evolving  wind  and  wave  conditions,  for  simple  fetch-limited  or 
duration-limited  wave  growth,  but  also  conditions  where  the  wind  direction  continuously  changes. 
This  generalization  of  TSA  can  be  considered  a  multiple-TSA  approach,  whereby  more  than  one 
broad-scale  parametric  term  is  allowed,  corresponding  to  multi  -peaked  spectra.  Thus,  this 
modified  TSA  can  handle  situations  where  the  wind  direction  veers,  or  in  swell-  windsea 
interactions  cases.  Model  accuracy  is  improved,  for  waves  generated  in  turning  winds,  or  by 
rapidly  developing  storms. 

b)  We  have  completed  a  variety  of  tests  and  model  implementations  with  the  new  generalized  TSA, 
involving  standard  modem  formulations  for  the  source  tenns  for  wind  input  energy,  and  wave 
dissipation  (Sin  and  Sds)  parameterizations,  including  state-of-the-art  so-called  ‘ST 4’ {Ardhuin  et 
al.,  2010).  Comparisons  were  made  with  more  historical  versions  of  these  source  terms,  for 
example  ‘ST1’,  dating  from  the  original  WAM  code.  We  have  shown  that,  without  any  tuning  of 
ST4  source  term  coefficients,  accurate  waves  simulations  are  possible  with  TSA  implemented  in 
WW3,  giving  the  best  overall  performance.  In  model  inter-comparisons,  implementations  were 
made  in  the  most  recent  versions  of  WW3  (4.10  and  4.18),  and  SWAN  (unstructured  grid,  and 
shallow  water  waves),  and  tests  were  conducted  to  simulate  waves  generated  in  stonn  conditions, 
again  showing  that  ST4  source  terms  and  TSA  can  demonstrate  the  best  overall  model  skill. 

c)  Optimized  versions  of  TSA  have  been  fonnulated,  in  terms  of  computer  mn  time,  efficiency,  and 
model  accuracy,  and  implemented  in  very  recent  versions  of  WW3  and  SWAN.  This  involves  an 
optimized  selection  of  grids  and  computation  loops  over  which  the  computation  is  carried  out, 

Thus  far,  for  implementations  in  WW3,  TSA  has  achieved  a  speed  that  is  approximately  20x 
slower  than  a  standard  DIA  fonnulation  for  the  nonlinear  wave-wave  interactions.  Tests  have 
included  waves  generated  by  hurricane  Juan  (2003),  as  well  as  several  mid-latitude  nor  'easters. 


IMPACT  /  APPLICATIONS 


The  value  of  this  project,  to  the  scientific  community  and  to  society,  is  that  a  more  physically  based 
wave  models  will  be  developed,  compared  to  the  present  state-of-the-art,  which  is  DIA,  for  the 
nonlinear  wave-wave  interactions,  at  the  heart  of  operational  wave  forecast  models.  TSA  has  been 
established  as  a  viable  (more  accurate)  alternative  to  DIA.  Thus,  where  inaccuracies  in  DIA  were 
hidden  in  tuning  coefficients,  for  example  embedded  in  Sin  and  5*,  it  is  now  be  possible  to  better 
simulate  and  model  the  processes,  previously  hidden  by  tuned  coefficients.  Wave  forecasts  should 
improve,  for  global,  basin,  regional  and  coastal  scale  studies,  with  application  of  TSA  as  the  key  part 
of  the  forecast  wave  model.  The  success  of  this  project  can  lead  the  way  for  TSA  being  implemented 
in  operational  wave  forecast  models,  such  as  WW3  and  SWAN,  defining  a  new  state-of-the-art. 

National  Security 

Landfalling  severe  storms  and  hurricanes,  like  Katrina  (2005)  or  Sandy  (2012)  generate  large  waves 
that  cover  the  coastal  and  nearshore  areas,  causing  significant  damage  and  sometimes  loss  of  life. 
Better  understanding  and  forecasting  of  these  storms  and  the  ocean  waves  that  they  can  generate,  can 
serve  to  potentially  reduce  some  of  these  problems.  This  project  focuses  on  the  task  of  obtaining  better 
science-based  parameterizations  for  ocean  wave  models,  suitable  for  operational  marine  forecasting. 

Economic  Development 

The  coastal  zone  involves  significant  potential  future  economic  development,  e.g.,  residences, 
recreation,  fisheries,  aquaculture,  coastal  transportation.  Better  forecasts,  with  longer  lead-time,  and 
better  accuracy  can  help  reduce  potential  risk  to  these  economic  developments,  due  to  ocean  waves. 

Quality  of  Life 

Development  of  the  coastal  zone  involves  residences,  recreation,  fisheries,  aquaculture,  coastal 
transportation.  Ocean  waves  are  a  risk  to  these  developments  and  activities,  both  from  storm  -related 
damage,  and  also  of  anthropogenic  factors,  such  as  transport  of  pollutants,  and  lastly,  future  potential 
issues  such  as  changing  wave  climate  and  wave  impacts  on  these  environments.  All  represent  potential 
future  impacts  on  Quality  of  Life,  e.g.,  public  and  ecosystem  health,  resource  management. 

TRANSITIONS 

The  new  parameterization,  TSA,  is  implemented  and  tested  in  WW3,  which  is  the  operational  forecast 
wave  model,  for  NOAA  /  NCEP  and  several  other  marine  forecast  offices,  internationally.  Thus,  when 
tests  are  completed,  which  is  scheduled  to  occur  imminently,  the  resulting  model  will  be  easily  useable 
as  the  operational  forecast  model  for  these  marine  forecast  offices. 

RELATED  PROJECTS 

A  related  project  is  supported  by  the  Marine  Environmental  Observation  Prediction  and  Response 
Network  (MEOPAR;  http://meopar.ca/)  a  team  of  outstanding,  Canadian  natural  and  social  scientists, 
working  to  better  understand  and  predict  the  impact  of  marine  hazards  on  human  activities  and 
ecosystems. .  .and  improve  response.  Perrie’s  MEOPAR  project  focuses  on  improved  wave  forecasts 
for  waters  off  the  coast  of  Atlantic  Canada,  with  particular  focus  for  Halifax  Harbor  and  approaches. 


REFERENCES 


Ardhuin,  F.,  E.  Rogers,  A.  V.  Babanin,  J.  F.  Filipot,  R.  Magne,  A.  Roland,  A.  van  der  Westhuysen,  P. 

Queffeulou,  J.  M.  Lefevre,  L.  Aouf,  and  F.  Collard  (2010),  Semiempirical  dissipation  source 

functions  for  ocean  waves.  Part  I:  Definition,  Calibration,  and  Validation,  J.  Phys.  Oceanogr., 

40(9),  1 9 1 7-1 94 1 ,  doi:  1 0. 1 1 75/20 1 0JPO4324. 1 . 

Perrie,  W.,  and  D.  Resio,  2009:  A  Two-Scale  Approximation  for  Efficient  Representation  of  Nonlinear 
Energy  Transfers  in  a  Wind  Wave  Spectrum.  Part  II:  Application  to  Observed  Wave  Spectra.  J. 
Physical  Oceanography.  Vol.  39,  pages  2451  -  2476.  DOI:  10.1 175/2009JPO3947.1 

Perrie,  W.,  Toulany,  B.,  Resio,  D.,  Roland,  A.,  and  J.-P.  Auclair,  2013:  A  two-scale  approximation  for 
wave-wave  interactions  in  an  operational  wave  model:  Ocean  Modelling.  70,  38-51. 
doi.org/10. 10 16/j.ocemod.20 13.06.008 

Resio,  D.,  Perrie,  W.,  2008:  two-scale  approximation  for  efficient  representation  of  nonlinear  energy 
transfers  in  a  wind  wave  spectrum.  Part  1:  Theoretical  Development.  J.  Phys.  Ocean.  38:  2801-2816. 

Tolman  H.  L.,  and  Chalikov  D.V.,  1996:  Source  terms  in  a  third-generation  wind-wave  model.  J.  Phys. 
Oceanogr.,  26,  2497-2518. 

WAMDI  Group  (1988).  The  WAM  model  -  a  third  generation  ocean  wave  prediction  model.  J.  Phys. 
Ocean.,  18,  1775-1809. 

PUBLICATIONS 

a.  Articles  in  refereed  publications 

1.  Chen,  Z.,  Zhang,  B.,  He,  Y.,  Qiu,  Z.,  Perrie,  W.,  2014:  A  new  Modulation  Transfer  Function 
for  the  Retrieval  of  Ocean  Wave  Spectra  from  X-band  marine  radar  imagery.  Submitted  to 
Chinese  Journal  of  Oceanology  and  Limnology. 

2.  Hwang,  P.,  Perrie,  W.,  and  Zhang,  B.,  2014:  Cross  Polarization  Radar  Backscattering  from  the 
Ocean  Surface  and  Its  Dependence  on  Wind  Velocity.  Geoscience  Remote  Sens.  Lett,  doi: 

10.1 109/LGRS. 2014.2324276 

3.  Li,  M.  Z.,  C.  G.  Hannah,  W.  Perrie,  C.  Tang,  R.  H.  Prescott  and  D.  A.  Greenberg,  2014: 
Modelling  Seabed  Shear  Stress,  Sediment  Mobility  and  Sediment  Transport  in  the  Bay  of 
Fundy.  Submitted  to  CUES  Canadian  Journal  of  Earth  Sciences 

4.  Liu,  G.  and  W.  Perrie,  2014:  Sea-State-Dependent  Wind  Work  on  the  Geostrophic  Ocean 
Circulation.  Geophys.  Res.Lett.,  40,  3150-3156,  doi:  10. 1002/grl. 50624. 

5.  Liu,  G.,  Perrie,  W.,  and  He,  Y.,  2014:  Ocean  surface  Stokes  drift  from  scatterometer 
observations,  International  Journal  of  Remote  Sensing,  35:5,  1966-1978,  DOI: 
10.1080/01431161.2014.880818 

6.  Perrie,  W.,  Toulany,  and  Resio,  D.,  2014:  Generalized  two-scale  approximation  for  nonlinear 
wave-wave  interactions  for  operational  wave  forecasting.  Almost  submitted  to  Ocean 
Modelling. 

7.  Shen,  H.,  W.  Perrie,  Y.  Hu,  S.  Li,  Y.  He,  2014:  Remote  sensing  of  waves  propagating  in  the 
marginal  ice  zone.  Submitted  to  J.  Geophys.  Res. 

8.  Sun,  Y,  Perrie,  W.,  Toulany,  B.,  2014:  Forecast  skill  for  the  two-scale  approximation 
implemented  in  unstructured  grid  SWAN  wave  model.  Almost  submitted  to  Weather  and 
Forecasting. 

9.  Xu,  F.,  Bui.  T.  T.  D.,  Perrie  W.,  2014  :  The  observed  analysis  on  the  wave  spectra  of  Hurricane 
Juan  (2003).  Accepted  by  Acta  Oceanologica  Sinica. 

10.  Xu,  F.,  X.  Xiao,  W.  Perrie,  2014:  Simulation  of  Cyclone-Generated  Waves  in  Coastal  Shallow 
Waters.  In  International  Society  of  Offshore  and  Polar  Engineers,  Journal',  7  pages. 


I 

1 0.5 

1™0.25 
0.125 
0.0625 
0.03125 
-  0.015625 
-  0.0078125 
-0.00390625 
-0.00195312 
-  0.000976562 
-  0.000488281 
-0.000244141 
-0.00012207 
-  6.10352e-05 


J3.05176*— 05 
1 .52588e-05 


n(ffy 


0.2 

f(Hz) 


(c)  2-d  action 


(d)  DIA 


(e)  WRT 


(f)  TSA 


0.5005 

0.25025 

1-0.125125 
0.0625625 
0.0312812 
0.0156406 
0.00782031 
0.00391016 
HO 


-0.00391016 

-0.00782031 

-0.0156406 

-0.0312812 

-0.0625625 

-0.125125 

-0.25025 


-0.5005 


Sni(f,Q) 


Figure  1.  (a)  Decomposition  in  terms  of  broad-scale  and  local-scale  terms  nonnalized  by  the  f~4,  (b)  1- 
d  comparison  of  DIA,  WRT  and  TSA,  (units:  m  ),  (c)  2-d  action  density  (d)  Sni(f,9)  from  DIA,  (e) 
WRT,  and  (f)  TSA.  Color-bar  for  DIA,  WRT,  TSA  scales  to  ±3. 36x1  O'5,  ±3.13xl0'6,  and  ±2.64x1  O'6, 
respectively.  Color-bar  for  n(f,9)  scales  to  x  10. 1.  The  1st  peak  is  at  0.0788  Hz,  the  2nd  at  0.135  Hz, 
shifted  90°  with  respect  to  1st  peak.  Both  distributions  are  JONSWAP:  y=3.3.  From  Perrie  et  al  (2013). 


